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FOREWORD 


This  interim  technical  report  was  prepared  by  the  School  of  Mechanical  Engin¬ 
eering,  Georgia  Institute  of  Technology,  Atlanta,  Georgia,  on  contract  AF  33(6l6)- 
8396  for  the  Aeronautical  Research  Laboratories,  Office  of  Aerospace  Research, 
United  States  AVr  Force.  The  report  contains  a  summary  of  the  results  obtained 
from  an  analytical  and  experimental  study  of  the  influence  of  resonant  acoustic 
vibrations  on  viscous  laminar  channel  flow  and  was  performed  on  Task  7063-01, 
"Research  in  Heat  Transfer  Phenomena"  of  Project  7063,  "Mechanics  of  Flight".  The 
research  was  administered  under  the  technical  cognizance  of  Dr.  Max  Scherberg  of 
the  Thermomechanics  Research  Laboratory,  ARL.  Dr.  Scherberg* s  interest  in  the 
work  and  his  many  helpful  suggestions  are  deeply  appreciated. 
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ABSTRACT 


This  report  presents  an  anauLytical  solution  for  low  Mach  nianber 
viscous  compressible  channel  flow  \inder  the  Influence  of  a  resonant  acoustic 
field.  The  solution  predicts  the  existence  of  standing  vortices  in  the  flew 
and  that  the  ratio  of  the  maximxm  width  of  the  vortices  to  the  half  width 
of  the  channel  is  solely  a  function  of  the  parameter  where  M  is  the 
Mach  number  based  on  the  average  through  flow  velocity  and  M^  is  the  acoustic 
Mach  number  based  on  the  maximm  amplitude  of  the  time  dependent  velocity. 

The  order  of  magnitude  analysis  which  was  used  to  reduce  the  equations  of 
motion  required  that  M/M^  shovild  be  less  than  approximately  ten. 

The  existence  of  stemding  vortices  was  confirmed  experimentally  by  a 
visualization  study  of  tube  flow.  Typical  flew  fields  are  shown  both 
graphically  and  photographically. 

Finally,  the  threshold  data  of  Eastwood,  ARL  62-326,  for  convective 
heat  transfer  londer  the  influence  of  a  resonant  acoustic  field  axe  correlated 
for  the  range  of  tube  diameter  Reynolds  numbers  between  11, OCX)  and  33 >000. 

It  was  shown  for  this  case  that  sound  affects  the  local  heat  transfer  rate 
when  M/M^  is  less  theui  63- 
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Average  kinematic  viscosity  of  a  fluid  ft  /sec 

p  Density  of  a  fluid  Ib^/ft^ 

Time  mean  component  of  density  Ib^/ft^ 
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CHAPTER  I 


INTRODUCTION 

A  program  to  study  the  effects  of  a  resonant  acoustic  field  on  viscous 
fluid  flow  and  convective  heat  transfer  was  initiated  in  the  School  of 
Mechanical  Engineering  of  the  Georgia  Institute  of  Technology  in  1958  under 
the  sponsorship  of  the  Aeronautical  Research  Laboratory  at  Wright  Field. 

Since  its  beginning,  the  study  [l,2]  has  shown  that  a  resonant  sound  field 
had  the  following  effects  on  convective  heat  transfer  in  a  4  inch  ID  horizontal 
isothermal  tube: 

1.  For  Reynolds  numbers  less  than  33»000  the  local  heat  transfer  co¬ 
efficient  varied  periodically  with  axial  position  in  the  tube. 

The  period  of  this  variation  was  one  half  the  wavelength  of  the 
impressed  resonant  acoustic  field.  Local  increases  and  decreases  of 
up  to  300  per  cent  occurred  in  the  heat  transfer  coefficient  with  the 
maximum  values  of  the  local  coefficient  occurring  at  the  sound  particle 
velocity  loops  and  the  minimum  values  occurring  at  the  velocity  nodes. 
Sound  threshold  values,  sound  pressure  levels  below  which  no  noticeable 
change  in  the  heat  transfer  characteristics  occurs,  were  found  to  in¬ 
crease  with  an  increase  in  Reynolds  number.  Only  slight  increases 

(30  per  cent)  in  the  overall  heat  transfer  coefficient  were  obtained. 
See  Figures  1  and  2  for  typical  results. 

2.  For  Reynolds  numbers  above  40,000  a  radical  change  in  the  local 
characteristics  took  place.  The  resonant  sound  waves  no  longer  in¬ 
creased  the  local  coefficients.  The  local  heat  transfer  coefficients 
still  varied  periodically  with  axial  position  and  the  period  was  the 
same  as  that  for  the  low  Reynolds  number  data.  The  position  of  the 
maximum  value  of  the  local  coefficient,  however,  shifted  to  the  node 
for  the  sound  particle  velocity  and  the  minimum  shifted  to  the  velocity 
loop.  This  shift  was  also  accompanied  by  a  radical  lowering  of  the 
sound  threshold  level.  The  overall  heat  transfer  coefficient  was  re¬ 
duced  in  all  cases  for  Reynolds  numbers  greater  than  40,000  and  less 
than  100,000.  See  Figure  3  for  typical  results. 

It  should  be  pointed  out  that  the  tube  for  which  the  above  results  were 
obtained  had  a  bell-mouth  inlet  so  that  the  flow  considered  was  not  fully 
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developed.  Calculations  indicate  that  the  boundary  layer  would  probably  be 
laminar  for  Reynolds  numbers,  calculated  using  the  diameter,  up  to  approximately 
33,000.  It  can  be  concluded  from  the  above  remarks  that  two  distinctly  different 
phenomena  exist;  one  characterized  by  laminar  flow  and  the  other  by  turbulent 
flow. 

In  order  to  achieve  a  better  understanding  of  the  phenomena  associated  with 
laminar  flow,  an  analytical  investigation  of  viscous  compressible  channel  flow 
under  the  influence  of  a  resonant  acoustic  field  was  initiated.  The  results  of 
this  analytical  study  are  presented  herein.  Preliminary  results  of  a  qualitative 
visualization  study  of  tube  flow  are  also  nresented  in  support  of  the  analytical 
interpretations. 
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CHAPTER  II 


ANALYTICAL  SOLUTION  FOR  CHANNEL  FLOW 

Aj _ Theoretical  Considerations 

The  problem  is  that  of  determining  analytically  the  velocity  field  which 
is  formed  by  imposing  a  resonant  acoustic  field  with  a  circular  frequency,  u, 
upon  a  compressible  fluid  flowing  between  infinite  parallel  walls  separated  by 
a  distance  2a.  The  mathematical  model  is  shown  in  Figure  4>  The  flow  is  assumed 
to  be  independent  of  the  z  direction.  Thus  the  velocity  of  the  fluid  is  depend¬ 
ent  upon  position  as  given  by  x  and  y  and  upon  time,  t.  For  this  model  the 
equations  which  govern  the  fluid  motion  are 
Continuity 


Pt  +  ~  ° 


...(1) 


Momentum 


[^(2  +  V^})]^  + 

*  ’'x’ly 

Y  -  p  + 

y 

C|>(2  Vj,  -  |(Uj,  +  Vy!)]y  + 

U(Uy  + 


...(2a) 


...(2b) 
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3.  Assumptions 


In  order  to  reduce  the  above  equations  to  a  form  which  can  be  solved,  the 
following  assumptions  are  made: 

1.  It  is  assumed  that  the  pressure  may  be  written  as  the  sum  of  a  time 
mean  component  p(x,y)  and  a  time  dependent  component  pi(x,t)  which  is 
independent  of  the  y-direction.  It  is  further  assumed  that  P2^(x,t)  is 
directly  rslated  to  the  resonant  acoustic  field  and,  therefore,  it  is 
periodic  in  both  time,  t,  and  position,  x.  Thus 

p(x,  y,  t)  =  p(x,  y)  +  p^(x,  t)  ...(3) 

2.  Since  the  pressure  is  assumed  to  consist  of  a  time  mean  component  and 
a  periodic  component  it  is  plausible  to  assume  that  the  density  may 
also  be  represented  in  this  fashion.  However,  spacial  variations  in 
the  density  are  assumed  to  have  only  a  small  effect  upon  the  flow  and, 
therefore,  the  time  mean  of  the  density,  P(x,  y),  is  considered  to  be 
constant  and  equal  to  p^.  Thus 

p(x,  y,  t)  =  +  p^(x,  t)  ...(4) 

3.  The  viscosity  of  the  fluid  is  assumed  to  be  constant  and  equal  to 

4>  Following  Lin  [3],  it  is  assumed  that  the  x-velocity  component, 

^(x,  y,  t),  and  the  y-velocity  component,  v(x,  y,  t),  may  be  separated 
into  time  mean  and  periodic  components  as  follows: 


u(x,  y,  t) 

=  u(x,  y)  +  u^(x,  y,  t) 

...(5) 

v(x,  y,  t) 

=  v(x,  y)  +  v^(x,  y,  t) 

...(6) 

where 

^1  =  ^1 

=  0 

... (7) 

5.  No  body  forces  are  assumed  to  be  present. 

i.e. 


X  =  Y  =  0 


•  •  • 


(&) 
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6.  The  continuity  and  momentum  equations  can  be  written  in  dimensionless 
form  by  assuming  the  following  change  of  variables! 


u' 

=  u/u^  =  u  +  u| 

.„.(9a) 

v' 

II 

> 

II 

<1 

+ 

< 

...(9b) 

x' 

=  xA 

...(9c) 

y' 

=  yA 

...(9d) 

p' 

=  p/p^uf  =  p'  +  Pi 

...(9e) 

P' 

=  p/pQ  =  1  +  Pi 

...(9f) 

p' 

o 

II 

...(9g) 

t' 

=  t/(\/u^) 

...(9h) 

7.  The  order  of  magnitude  of  each  of  the  dimensionless  quantities  is 
assumed  to  be 


x* 

= 

e(i) 

....(lOa) 

y' 

= 

0(6) 

...(10b) 

4 

= 

e(i) 

...(10c) 

U 

= 

0(6) 

...(lOd) 

''i 

= 

0(6) 

...(lOe) 

V 

= 

0(6^) 

...(lOf) 

Pi 

= 

9(5) 

...(lOg) 
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C.  Reduced  Equations 


Continuity 

The  dimensionless  continuity  equation  can  be  reduced  by  retaining  only 
the  terms  which  are  of  order  one.  This  yields  the  following  dimensional  equa¬ 
tion 


p,  +  p  .u,  +  p  V,  =  0 

*^1^  ^0  1^  *^0  1. 


...(11) 


Taking  the  time  average  over  one  period  yields 


u  +  V  = 
X  y 


-p  [Pr“i  +  Pi-^i  +  } 


...(12) 


Momentum 

The  dimensionless  momentum  equation  can  be  reduced  by 

(1)  eliminating  the  pressure  between  the  x  and  y  equation  by  noting  that 

P  =  P  f 
'^xy  "^yx’ 

(2)  assuming  the  Reynolds  number,  PqUq^/Pq»  Is  of  order  lA  ,  snd 

(3)  retaining  only  the  terms  that  are  of  order  lA  or  larger,  or 

(4)  taking  the  time  average  and  retaining  only  the  terms  that  are  of 
order  one  or  greater. 

Conditions  (1),  (2)  and  (3)  reduce  the  momentum  equations  to 


"o  '^l 


=  tu,  +  Ui-Ui  +  V  .U  ] 

yyy  t  X  y 


Conditions  (1),  (2)  and  (4)  reduce  the  momentum  equations  to 


...(13) 


V  u 

0  yyy 


where 
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...(15) 


Q 


f  ^  ^  ^ 

■<U,.U,  +V,  .u,  +  “*-.U,  +-*’U,.U.  +  -* 

1.11  11  pi.  oil  p 

X  y  ^0  t  *^0  X  ^0 


Summary  of  the  Reduced  Equations 

(1)  For  u^(x,  y,  t)  and  v^(x,  y,  t) 

Continuity  ''i  “  ^^/Po^Pl 

■‘^x  y  '''t 


Momentum 


"o  '^l 


yyy 


+  U, 


+ 


...(16) 

...(17) 


The  solution  to  the  above  equations  must  satisfy  the  following  boundary 
conditions! 


Uj^  =  =  0  at  y  =  0  ...(l8a) 

Uj^  =  0,  =  0  at  y  =  a  ...(l8b) 

y 

Uj^  =  -U^  cos(ut)>  at  X  =  0,  y  =  a  ...(18c) 

The  first  and  second  conditions  are  required  by  no-slip  and  symmetry 
respectively.  The  third  condition  is  imposed  by  the  assumed  acoustic  field. 
Sanders  [4]  has  shown  experimentally  that  outside  a  very  thin  layer  of  fluid 
near  the  walls  the  velocity  is  the  same,  at  least  to  a  first  approximation,  as 
that  of  an  inviscid  compressible  fluid  undergoing  resonant  acoustic  vibrations. 

(2)  For  u(x.  v)  and  v(x.  v) 

Continuity  ^  +  v  =  -l/p  fp^^.u.  +  +  u^.p^  ]  ...(19) 

•  *  V  \/  V 


Momentum  v  u  =  Q  ...(20) 

0  yyy  y 

The  solution  to  the  above  equations  must  satisfy  the  following  boundary 
conditions : 
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u 

=  V 

s 

0  at  y  =  0 

...(21a) 

u 

y 

=  0, 

V 

=  0  at  y  =  a 

...(21b) 

f  Po^ 

^0 

aU 

=  constant  for  all  x 

...(2ic) 

The  first  two  conditions  are  again  no-slip  and  symmetry  respectively.  The 
third  condition  requires  steady  time  mean  flow  and  forms  the  basic  difference 
between  this  study  and  the  previous  studies  of  Rayleigh  [5]  and  Westervelt  [6] 
who  treated  the  special  case  of  no  through  flow. 

The  solutions  of  these  equations  will  be  outlined  in  detail  since  the 
underlying  assumptions  are  not  clearly  defined  in  the  existing  literature. 

D.  Solution  for  the  Time  Dependent  Velocity 

A  solution  to  Equations  I6,  17  and  18  can  be  obtained  by  the  method  of 
successive  approximations  [7],  i.e.,  the  velocities  u^(x,  y,  t)  and  v^(x,  y,  t) 
are  assumed  to  consist  of  two  terms  such  that 

Uj(x,  y,  t)  =  Ujq(x,  y,  t)  +  u^^Cx,  y,  t)  ...(22) 

Vj(x,  y,  t)  =  Vjq(x,  y,  t)  +  Vjj(x,  y,  t)  ...(23) 

where  the  first  approximations,  u^^  and  Vj^q  are  the  solution  of 

Continuity  u^q  ~  •'•(24) 

X  y  ■‘^t 

Momentum  v  u^q  =  Uj^q  +  F(x,  t)  ...(25) 

°  yy  t 

Equation  25  is  the  result  of  integrating  Equation  17  with  respect  to  y.  It 
can  also  be  obtained  from  an  order  of  magnitude  analysis  of  the  dimensionless 
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form  of  Equation  2a.  This  analysis  shows  that  F(x,  t)  = 


r  Pi ' 

0  X 


The  second  approximations,  and  are  the  solution  of 

Continuity  u.^  +  =  0 

^  X  y 


Momentum  u, ,  -  v  u, , 

11.  o  11 
t  yy 


~  "'^10‘'^10  "  ''io‘'^io 

X  V 


^(26) 

(27) 


The  solutions  to  the  above  equations  must  satisfy  the  following  boundary 
conditions! 


'^lO  “ll  “  ''lO  ‘  ''ll 

=  0 

at 

y 

=  0 

,...(28a) 

^0  =  '^ll  =  ''lO  =  ''ll 

=  0 

at 

y 

=  a 

...(28b) 

y  y 

u^Q  =  -U^  cos(ut),  = 

0  at 

y 

= 

a,  X  =  0 

„ ,.(28c) 

First  Approximations 

In  addition  to  the  seven  assumptions  given  previously  the  following 
assumptions  are  madei 

8.  It  is  assumed  that  all  quantities  which  are  periodic  in  time  are 
proportional  to  exp(-iut),  i.e 

f(x,  y,  t)  =  fj^(x,  y) 'expC-inwt)  ...(29) 

where  n  =  1,  2,  3,...  ..(30) 

9.  The  time  dependent  component  of  pressure,  pi(x,  t),  is  assumed  to  be 
the  same  as  the  pressure  in  a  stationary  inviscid  fluid  which  is  under¬ 
going  resonant  acoustic  vibrations.  Morse  [8]  gives  this  as 

Pj^(x,  t)  -  -i  c  sin{ux/c)  exp(-iu)t)  ...(31) 

10.  It  is  further  assumed  that  the  form  of  the  first  approximation  to 
the  x-velocity,  Uj^,  is 
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Ujq(x,  y,  t) 


...(32) 


=  Uq  cos(*~)  Fj(y)  exp(-iut) 

This  form  is  selected  since  Sanders  [4]  has  shown  that  a  real  fluid 
undergoing  resonant  acoustic  vibrations  has  the  same  velocity  as  an 
inviscid  fluid  outside  a  very  thin  layer  of  fluid  near  the  stationary 
boundary.  The  velocity  for  an  inviscid  fluid,  as  given  by  Morse[8], 
is 


u.(x,  t)  =  -U  cos(^)  exp(-iu)t) 
X  o  c 


...(33) 


Thus  the  fact  that  the  fluid  is  viscous  is  assumed  to  affect  only 
the  amplitude  and  the  phase  of  the  velocity  and  to  be  solely  a 
function  of  y,  namely  Fj(y). 

11.  Finally,  in  order  to  eliminate  p]^(x,  t)  from  the  above  equations,  it 
is  necessary  to  assume  a  relationship  between  the  pressure  and  the 
density  of  the  fluid.  The  pressure^ensity  variation  due  to  the  time 
mean  flow,  i.e.  due  to  p(x,  y)  and  f(x,  y),  will  be  omitted  from  this 
analysis  and  only  the  relationship  between  pressure  and  density  during 
the  compressions  and  rarefactions  of  the  fluid  will  be  considered.  To 
achieve  this  it  is  assumed  that 

(a)  The  fluid  is  a  perfect  gas,  i.e. 


p  =  Cp  *  Constant}  and 
(b)  The  compression-rarefaction  process  is  isentropic,  i.e. 
pp  =  Constant,  where  Ic  =  Cp/c^. 


The  wave  propagation  speed,  c,  for  isentropic  compressions  and 
rarefactions  is  given  by 


c2  =  2il  I  ^  ...(34) 

3p  's  p 

From  these  assumptions  and  the  further  assumption  that  p^/p  =^*9(&), 
the  following  relationship  can  be  shown 


Pi 


...(35) 
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Under  the  above  assumptions  Equation  25  reduces  to 


[fJI  +  iu/v^  +  iu/v^]  U^  cos(iJx/c)  exp(-i(Jt)  =  0  ...(36) 

For  this  to  be  satisfied  in  general  the  equation  in  brackets  must  be  equal  to 
zero,  i.e. 

F"  +  iu/v  F  =  -  Wv„  ...(37) 

1  0  1  0 

The  general  solution  for  Fj^(y)  is 

Fj^(y)  =  -1  +  A  exp[p(l-i)y]  +  B  exp[-p (l-i)y]  ...(38) 

where  p  =  ^u/2*v^ 
or  finally 

Ujq(x,  y,  t)  =  cos(ux/c)[-l  +  A  exp[p(l-i)y]  +  B  exp[-p(l-i)y]]  ...(39) 

The  y  component  of  velocity,  v^^q,  can  now  be  determined  from  Equation  24  since 
Pj^(x,  t)  and  'Jj^q(x,  y,  t)  are  known.  Upon  applying  the  appropriate  boundary 
conditions  to  the  real  parts  of  u^^  and  v^^q,  the  following  solution  is  obtainedt 

"  "^o  cos(ux/c)[cos(wt)  - 

®xp(-py)*cos(ut  -py)]  ...(40) 

U  u 

v,n(x,  y,  t)  =  sin(ux/c)[cos(ut  -  it/4)  - 

exp(-^y)*cos(ut  -py  -  n/4)]  ...(41) 

'The  first  approximations  to  u^  and  v^  given  above  are  the  same  as  those 
given  by  Westervelt  [6]. 
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Second  Aporoxlmations 


The  following  equation  is  the  result  of  substituting  the  complex 
form  of  Equations  40  and  41  into  Equation  27t 

Ujj  -  u^j  =  ^  sin{2ux/c).[l  -  exp(-p[l-i]y)]  exp(-2iut)  .«o(42) 
t  yy 

Another  assumption  is  now  necessary. 

12.  It  is  assumed  that  the  form  of  given  by 

uU^ 

y,  t)  =  sin(2ujx/c)  F2(y) •exp(-i2cjt)  ....(43) 


The  above  assumption  is  based  upon  the  premise  that  the  convective  terms 
give  rise  to  higher  harmonics  in  both  position  and  time. 

Substituting  Equation  43  into  Equation  42  produces  a  linear  total  dif¬ 
ferential  equation  in  of  Equations  26  and  42  under  the 

restraints  of  Equations  28  is 

y»  t)  ~  sin(2ux/c){£in{2ut)  + 

exp(- y2  py)*sin(2ut  -  >/2  py)  - 

2  exp(-  py)*sin(2u)t  -  py)}  .„.i44) 


and 


Vji(x,  y,  t) 


uU^ 

-  — ^  cos(2ux/c)[py  sin(2u)t)  + 

2^0 

(^2  -  1/2)  cos(2ut  +  it/4)  +  1/2  exp(- V2  py)- 

cos(2ut  -  ^2  py  +  n/4)  -  ^2  exp(-  py)" 

cos(2u)t  -  py  +  x/4)]  .0.(45) 
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Dimensionless  Velocity  Components 


The  following  two  equations  represent  the  time  dependent  components 
of  velocity,  Uj^  and  Vj^,  in  dimensionless  form  to  a  second  approximation? 

u^(x',  y',  t')  =  -  cos(2itx ' ){  cos(2nt' /M^)  - 

exp(- y')  cos(2itt'/M^  ” ^ 

(M^A)  •sin(4^x'  )[sin(4’tt'/Mj^)  + 
exp(-  A/2Rg  y' )  •sin(4’tt' /Mq  ~ ' 

2  exp(- ^R  y’ )  ■■sin(4)it'/M„  -  v^R^  y')j  ...  .(4'') 

and 

v^(x',  y',  t')  =  (-1/2  n/  yR^)  sin(2itx')'’ 

fcos(2itt'/M  -  itA)  -exp(-yR  y' )  “Cos(2itf/M^  -yP.  y  -itA)( 
-  (it/  yRg)  cos(4itx')[  yRg  y'  sin(4’tt'/MQ)  + 

(  y2  -  1/2)  cos(4itt'/MQ  +  it/4)  +  1/2  exp(-  y^F  y  )■ 
cos(4itt'/M^  " exp(-^/Rg  y')' 
cos(4itt'/M|^  "  ...(47) 

where 


R 


s 


=  21i4 

V 

o 


and 


(48) 


M 


o 


. . ..  (49) 


The  above  solution  is  valid  only  if  the  original  order  of  magnitude 
assumptions  are  satisfied.  Therefore,  the  following  restrictions  are  required: 
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...(50) 


(a)  =  9(6)  or  =  9(6) 

(b)  Rg  =  9(1/6^)  ...,(51) 


The  first  requirement  was  also  necessary  in  the  development  of  Equation  35. 

The  second  requirement  is  satisfied  by  gases  at  normal  pressures  and  tempera¬ 
tures  and  for  wavelengths  which  are  greater  than  approximately  6  inches. 

It  can  also  be  seen  that  the  contribution  of  the  second  approximation  is 
directly  proportional  to  the  magnitude  of  the  pressure  disturbance. 

Under  the  order  of  magnitude  assumptions  and  for  the  range  of  T^equency 
that  was  investigated  experimentally  (600  to  l600  cps),  it  can  be  shown,  numer¬ 
ically  that  the  second  approximation  can  be  neglected.  Therefore,  the  analysis 
of  the  time  mean  flow  will  be  based  on  the  first  approximation  to  the  time 
dependent  velocity,  u^q.  Uj^(0,  y',  t')  is  shown  in  Figure  5. 

E.  Time  Mean  Velocities 


The  primary  characteristics  of  the  time  mean  flow  can  be  determined  by 

considering  only  the  first  approximation  to  u^  and  v^,  i.e.  Equations  40  and 

41.  A  further  assumption  is,  however,  necessary. 

Pi  Pi 

13.  The  terms  and  v. 'U^  will  be  neglected  since  they  yield 

Po  Po  V  y 

terms  which  are  higher  harmonics  in  both  time  and  position. 
Therefore,  to  a  first  approximation 


0  =  “10-“10 

X  Y  ^0  t 


.(52) 


or  upon  substituting  Equations  40  and  41  into  Equation  52 


uU^ 

Q  =  -  •jf'  sin(2u)x/c)[2  -  exp(-Py)«[3  cos(Py)  - 
sin(Py)]  +  exp(-2Py)} 


(53) 
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Equation  20  can  now  be  solved  for  u(x,  y)  as  follows: 


Substitution  of  Equation  53  into  Equation  20  yields 

Vy  "  ^  sin(2Dx/c)[^  [exp(-Py)- 

(3  cos{Py]  -  sin{3y])  -  exp(-?Py)]]  ...(54) 

Equation  54  is  a  linear  partial  differential  equation  of  the  form 

Uyyy  =  f(x)*g(y)  ...(55) 

y  y'y" 

f(x)  I  I  ^  g(y)dy'dy"dy"'  + 

2 

F3(x)^  +  F^(x)y  +  F^(x)  ...(56) 

^  2 

sin(2tjx/c)(-4(py)^  +  2  ^y  +  3  - 

exp(-Py)[2  cos(py)  +  6  sin(Py)  +  exp(-Py)]}  + 

F^{x)‘y^/2  +  F^(x)*y  +  Fj(x)  ...(57) 

From  the  no-slip  condition,  F^(x)  =  0,  the  y  component  of  velocity,  v, 

can  now  be  determined  from  continuity.  Equation  19.  This  yields 

U?u  2 

V  =”““2  cos(2ux/c)[-4(Py)‘^  +  2py  +  3  - 
^  4c 

exp(-Py)[2  cos(Py)  +  7  sin(Py)]  -  exp(-^y)}  + 

U^u  ,  2  , 

exp(-Py)  sin(py)  -  Fj(x)  ^  -  F^(x)  y  ...(58) 


for  which  the  solution  is 
j(x,  y)  = 


or 


u(x,  y)  = 
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which,  w^en  integrated  with  respect  to  y  and  subjected  to  the  no-slip  condition, 
becomes 


v(x,  y)  =  - 


{[-§  (Py)^  +  2(py)2  +  6py 


10  + 


exp(-Py)[9  cos(Py)+5  sin(Py)} +exp(-2py)]  cos(2ux/c)  + 
exp(-Py)[cos(Py)  +  sinOy)}  -  ij-  -  Fj(x)  yV6-F^(x)yV2 


...(59) 


If  F^(x)  and  F^(x)  are  determined  from  the  boundary  conditions  given  in 
Equations  21b  and  21c,  then  the  time  mean  velocities  in  dimensionless  form  become 


u'(x',  y')  =  3(u/u^)Ly  -  Y^/2]  -  (M^S)  sin(4itx')* 

{[9  -  ^][y  -  yV2]  -  3  +  exp(-a)* 

[2  cos (a)  +  6  sin(a)  +  exp(-a)]]  ...(60) 

and 

V  (x',  y')  =  -  (M^ii/4)(aA){cos(4Jtx')[{9  -  -  Y^}  + 

6y  -  lO/pa  +  (l/pa)  exp(-a)‘[9  cos(a)  +  5  sin(a)  +  exp(-a)}]  + 
(l/pa)  exp(-a){cos(a)  +  sin(a)  -  l/pa]  ...(61) 


where  y  =  y/®  a  =  Py 


For  gases  at  normal  pressures  and  temperatures  p»l  and,  therefore,  l/pa«l 
and  can  be  neglected. 

The  above  solution  further  clarifies  the  order  of  magnitude  assumptions 
since  if  the  above  solution  is  valid  then 


(a) 

U/U^  =  9(6) 

or  M/M^  =  9(6) 

...(62) 

(b) 

=  0(5) 

which  from  (a)  yields 

...(63) 

M  =  9(5^) 

...(64) 
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and 


(c)  aA  =  e(6)  ...(65) 

The  X  component  of  the  time  mean  velocity  profile  is  shown  in  Figure  6  for 
selected  values  of  x'. 

£i _ Stream  Function 

If  the  special  variations  in  the  time  mean  density  are  neglected  then  an 
expression  for  the  stream  function  can  be  determined  as  follows t 

y* 

y*)  =  tj)(x',  0)  +  f  p  u(x',  z)dz  ...(66) 

o  ° 

If  tj)(x',  0)  =  0  then  \|)(x',  y)  becomes 

f(x',  y»)  =  4(x',  v^)/p^aU  =  |(y^  ~  yV3)  +  (mJ/m)  sin(4itx')- 

I”  2  “  yV3]  +  3y  +  (l/^a)[exp(“a)  • 

[8  cos(a)  t  4  sin(a)  +  exp(-a)}  -  9]]  ...(67) 

Figure  7  is  a  graphical  representation  of  Equation  67  for  selected  values  of 
3 1  X  9  M  snd  c 

If  only  one  half  of  the  flow  shown  in  Figure  7  is  considered  then  there  are 
three  vortices  for  each  increment  of  length,  \/2.  Two  of  the  vortices  are  very 
thin  and  would  be  undetected  in  flow  visualization  studies.  However,  the  third 
vortex  can  completely  fill  the  half  width  of  the  channel  if  the  through  flow  is 
reduced  to  zero.  This  condition  would  also  produce  a  fourth  vortex.  An  approxi¬ 
mate  expression  can  be  determined  for  the  size,  of  this  vortex  in  its  visible 
range  since  exp(-a)  0  and  also  l/^a  0.  The  size,  Yg*  is  simply  the  largest 
positive  root  of  i|)(x',  y*)  =  0.  Under  the  above  simplifications  Yg  becomes 
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Figures  8  and  9  ara  the  graf^ical  representations  of  Equations  70  and  71 
respectively. 
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CHAPTER  III 


FLOW  VISUALIZATION  STUDY 


A.  GeoiMtrv 

Before  the  analytical  study  presented  in  Chapter  II  was  initiated  a 
preliminary  flow  visualization  study  was  conducted  to  determine  vdiether  or  not 
a  resonant  acoustic  field  would  produce  a  time-mean  secondary  flow.  Since  the 
heat  transfer  experiments  had  been  conducted  for  a  circular  tube  it  seemed 
logical  to  use  a  circular  tube  for  the  preliminary  flow  visualization  study. 

The  first  results  of  this  study  (these  are  qualitative  only)  demonstrated  that 
secondary  time-mean-flow  does  exist  and  that  the  flow  seemed  to  be  periodic  with 
respect  to  axial  position  with  a  period  equal  to  one  half  the  wavelength  of  the 
impressed  standing  wave  field.  The  flow  patterns  are  shown  in  Figure  10. 

This  study  also  showed  that  when  secondary  flows  were  present  the  average 
through-flow  velocity  was  small  in  comparison  to  the  velocity  of  a  fluid  particle 
at  a  velocity  loop  in  the  standing  wave  field.  It  also  indicated  that  the 
pressure  disturbance  was  small  in  comparison  to  the  time-mean  pressure.  These 
observatlonsi  the  analyses  of  Rayleigh  [5]»  Westervelt  [6J,  and  Schlichting  [7] 
and  the  flow  visualization  studies  of  Andrade  [9]  formed  the  basis  of  the  analy¬ 
sis  which  are  presented  in  Chapter  II. 

It  was  originally  anticipated  that  a  study  of  the  two-dimensional  channel 
would  produce  the  simplest  experimental  system.  Consequently,  the  theoretical 
work  was  based  upon  the  infinite  parallel  plate  geometry.  This  geometry  also 
afforded  a  check  on  Rayleigh's  solution  for  the  special  case  of  no  through  flow. 

A  preliminary  investigation  of  a  rectangular  duct  verified  the  premise  that 
standing  vortices  could  also  be  produced  in  this  flow.  A  typical  vortex  is 
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shown  in  Figure  11.  However,  this  geometry  from  an  experimental  aspect  has  many 
iisadvantages.  Several  of  these  are  discussed  below. 

1.  In  order  to  improve  a  flow  visualization  apparatus  fiom  one  which  simply 
verifies  the  existence  of  standing  vortices  to  one  which  will  yield 
quantitative  data,  a  great  deal  of  care  must  be  taken  in  the  design  to 
be  certain  that  the  actual  flow  is  representative  of  the  postulated 
channel  flow.  Since  a  side  boundary  layer  is  formed  in  a  rectangular 
duct,  the  duct  width  to  height  ratio  should  be  quite  large  so  that  the 
side  boundary  layer  effects  can  be  neglected. 

2.  In  order  to  verify  the  analytical  solution,  the  easiest  parameter  to 

measure  is  the  size  of  the  vortices.  This  should  be  solely  a  function 
of  Accurate  measurements  from  photographs  require  a  duct  height 

of  at  least  0.5  inch.  Consequently,  a  duct  width  of  at  least  5  to  10 
inches  is  required. 

3.  In  order  to  supply  the  large  quantities  of  acoustic  power  required  with 
a  simple  driver-horn  system,  the  horn  should  also  be  at  least  0.5  inch 
in  diameter.  This  in  turn  limits  the  minimum  frequency  to  approxi¬ 
mately  1,000  cps  and  yields  a  sound  wavelength  for  an  open  duct  in 
normal  atmospheric  air  to  about  14  inches.  Thus  the  duct  width  is  of 
the  same  order  of  magnitude  as  the  sound  v;avelength. 

4.  In  order  to  have  plane  wave  propagation,  it  is  necessary  that  the 
largest  cross  sectional  duct  dimension  be  small  in  comparison  to  the 
wavelength.  Thus  a  contradiction  in  requirements  has  resulted.  The 
easiest  requirement  to  relax  is  the  duct  width.  However,  it  is  then 
questionable  whether  or  not  the  data  will  be  of  quantitative  value. 

5.  In  order  to  overcome  the  above  difficulties  a  new  method  of  producing 
the  sound  fiel.  would  have  to  be  devised. 

If  the  same  factors  are  considered  for  a  circular  tube  then  it  is  found 
that  no  outstanding  difficulties  exist  other  than  relating  the  quantitative 
results  for  a  tube  to  the  analytical  solution  for  a  channel.  Circular  tube 
flow  is  axially  symmetric  (no  extraneous  side  boundary  layer  effects).  Also, 
the  tube  diameter  can  be  made  small  in  comparison  to  the  sound  wavelength  with¬ 
out  sacrificing  accuracy  in  the  vortex  measurements. 

On  the  premise  that  the  vortex  size  should  be  a  function  of  the  same 
parameter  M/M^  for  both  a  circular  tube  and  channel  flow,  it  was  decided  to 
attempt  to  perfect  the  circular  tube  geometry  in  order  to  obtain  quantitative 
flow  visualization  data. 
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B.  Apparatus 

The  experimental  apparatus  consisted  of  the  following  items: 

1.  Test  Section.  The  test  section  was  composed  of  an  inlet  plenum 
chamber  and  a  1  inch  ID  Pyrex  glass  tube  42  inches  long.  For  the 
air  pressure  and  temperature  that  was  maintained,  the  fundamental 
resonant  frequency  for  this  tube  was  160  cps.  The  inlet  end  of  the 
tube  was  flared  to  form  a  bellmouthed  entrance.  Primary  air  was 
admitted  to  the  inlet  plenum  from  a  small  calming  chamber.  The  air 
flow  was  directed  up  and  to  the  rear  of  the  inlet  plenum.  This 
permitted  sufficient  residence  time  for  the  air  to  reach  thermal 
equilibrium  with  the  system  and  also  prevented  appreciable  distortion 
in  the  velocity  profile  at  the  tube  inlet. 

Smoke  was  also  admitted  to  the  inlet  plenum.  However,  the  smoke 
was  added  periodically  to  this  chamber  and  allowed  to  mix  with  the 
steady  flow  of  primary  air.  The  smoke  remained  in  the  inlet  plenum 
long  enough  to  be  completely  mixed  with  the  air  and  a  charge  of 
smoke  often  lasted  for  more  than  20  minutes.  The  air-smoke  mixture 
was  discharged  from  the  tube  into  an  exit  plenum  which  was  provided 
with  an  exhaust  blower  and  duct.  The  test  section  was  housed  in  an 
acoustic  shield  to  reduce  the  sound  level  in  the  laboratory  and  also 
to  eliminate  extraneous  light.  The  entire  system  is  shown  schematically 
in  Figure  12  and  photographically  in  Figure  13. 

2.  Light  Source.  Two  light  sources  were  used  in  this  study;  a  400  watt 
mercury  vapor  lamp  and  a  flash  tube  which  produced  a  200  watt-second 
line  source.  The  flow  was  illisninated  py  directing  the  light  beam 
down  the  tube  from  the  inlet  end.  Either  two  collimating  slits  were 
used  to  produce  a  light  beam  to  illuminate  a  vertical  plane  of  fluid 
through  the  axis  of  the  tube  or  the  beam  was  allowed  to  illuminate 
the  entire  tube. 

3.  Air  and  Smoke  Supply.  A  variable  speed  Roots  blower  supplied  low 
pressure  atmospheric  air  to  the  system.  The  aix  was  taken  directly 
from  the  air  conditioned  laboratory.  The  thermal  ineztia  of  the  system 
was  sufficient  to  damp  small  oscillations  in  the  room  temperature. 

Since  a  positive  displacement  blower  was  used  it  was  necessary  to  in¬ 
sert  a  storage  tank  in  the  supply  line  to  damp  the  low  frequency  flow 
pulsations.  The  air  flow  from  the  storage  was  measured  with  a  dry  gas 
meter  and  it  was  then  split  into  a  primary  air  stream  and  a  secondary 
stream  to  supply  air  to  either  a  tcbacco  smoker  or  to  an  oil  vaporizer. 
The  two  streams  were  united  again  in  the  inlet  plenum  chamber.  The 
flow  patterns  were  qualitatively  the  same  for  both  the  oil  and  the 
tobacco  smoke.  The  oil  unit  was  much  less  trouble  to  use  and,  there¬ 
fore,  it  constituted  the  primary  source  of  smoke  for  the  system. 

Large  variations  in  the  ratio  of  smoke  to  air  did  not  appear  to  affect 
the  flow  patterns,  unless  extremely  large  amounts  of  smoke  were  used. 

4.  Sound  Generating  Eouinment.  A  Hewlett-Packard  Model  206a  audio  signal 
generator  provided  a  source  of  continuously  variable  audio  frequency 


21 


voltage  at  a  total  distortion  level  of  less  than  0.1  per  cent.  This 
was  amplified  by  a  Bogen  75  watt  Challenger  amplifier.  The  amplifier 
output  was  measured  by  a  John  Fluke  Model  102  VAW  meter  and  was  then 
introduced  into  the  test  section  by  either  a  conical  or  a  cylindrical 
horn  coupled  to  a  University  type  PA-HF  speaker  driver  wb.ich  had  a 
frequency  response  of  70  to  10,000  cps .  A  wave  analysis  of  the 
resonant  acoustic  field  showed  a  distortion  level  of  less  than  one 
per  cent.  In  order  to  produce  an  axially  symmetric  flow  field  it  was 
necessary  to  maintain  a  very  low  distortion  level  in  the  tube.  The 
driver-horn  assembly  was  rack  mounted  so  that  the  position  of  the 
horn  in  the  tube  could  be  easily  varied.  This  adjustment  was  very 
critical,  especially  for  the  higher  frequencies  (960  to  1600  cps),  A 
schematic  diagram  of  the  sound  generating  system,  is  shown  in  Figure  12. 

5.  Sound  .Measuring  Eauipment ■  The  sound  measuring  equipment  included  the 
following  General  Radio  items t 


_ Ilsm _  _ lYPg 

Sound  level  meter  1551-A 

Power  supply  1262-A 

High  level  microphone  assembly  1551-PlH 

20  db  attenuator  pad  1551-P  11 

Sound  level  calibrator  i552-B 

Transistor  oscillator  1307-a 


In  addition  to  the  General  Radio  items  a  Hewlett-Packard  Model 
302-A  wave  analyzer  with  a  Hewlett-Packard  \todel  297-a  sweep  drive  and 
a  Mosley  Model  2-D  X-Y  recorder  were  used  to  analyze  and  record  the 
pressure  in  the  test  section. 

In  order  to  determine  the  axial  sound  pressure  districution  in 
the  tube  the  high  level  microphone  assembly  was  equipped  with  a  probe 
a.nd  was  mounted  on  a  trac.k.  A  ten  turn  potentiometer  and  a  dial  cord 
system  were  used  as  a  direct  current  voltage  input  to  tne  X  coordinate 
of  the  recorder,  This  voltage  varied  linearly  with  the  axial  position 
of  the  probe.  The  output  from  the  sound  level  meter  was  the  input  t.: 
the  Y  coordinate.  Ma.nual  movement  of  tne  pressure  probe  resulted  in 
a  continuous  plot  of  sou.nd  pressure  level  versus  axial  position., 

Figure  14  shows  a  typical  pressure  distribution,  /'ia-r  analyses  of  the 
resonant  acoustic  field  showed  that  the  distortion  level  was  less  than 
one  per  cent. 

A  Sorensen  Model  2501  A-C  line  voltage  regulator  supplied  the 
power  for  all  of  the  electronic  equipment, 

c.  Photographic  Equipment,  A  4-  x  5-inch  speed  press  ca.mera  with  a 

Schneider  Xenotar  1:2,  3,150  lens  a,nd  a  Pclaroid  Land  Oa.Tera  back  and 
30C0  speed  Polaroid  film  were  used  to  obtain  the  photographic  data. 
Both  time  and  flash  exposures  were  used. 
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A  warm  up  time  of  at  least  3  hours  was  provided  for  all  of  the  equipment. 

After  the  warm  up  period  the  following  procedure  was  followed: 

1.  The  iesired  flow  rate  was  set  up  using  only  the  primary  air  supply 
(no  flow  through  the  smoke  generator). 

2.  The  desired  resonant  frequency  was  set  on  the  audio  signal  generator. 

3.  The  amplifier  gain  was  increased  until  the  desired  sound  level  was 
reached.  The  maximum  sound  pressure  level  was  determined  by  insert¬ 
ing  the  calibrated  probe  into  the  test  section, 

4. -  The  driver-horn  assembly  was  then  adjusted  to  minimize  the  driver 

power.  This  was  a  very  accurate  and  simple  method  of  tuning  the 
system.  After  this,  additional  fine  adjustments  in  the  frequency  and 
the  horn  position  were  made  until  the  power  reached  a  relative  minimum, 

5.  The  probe  was  then  removed  after  the  maximum  sound  pressure  level 
was  recorded. 

6.  Photographs  of  the  flow  field  were  taken  after  allowing  sufficient 
time  for  steady  state  conditions  to  exist. 

7.  The  data  that  were  recorded  are:  frequency;  power  supplied  to  the 
driver;  maximum  sound  pressure  level  and  flow  rate,  temperature  and 
pressure  of  the  air. 

From  the  data  obtained  the  vortex  size  and  the  value  of  the  parameter 

^  2 

m/m  were  determined. 

0 
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CHAPTER  IV 


DISCUSSION  OF  RESULTS 

At  the  initiation  of  this  study  it  was  thought  that  an  experimental  flat 
plate  geometry  would  give  simpler  flow  and  heat  transfer  characteristics. 
Consequently,  an  analytical  and  experimental  study  was  initiated  utilizing  this 
geometry  and  its  relation  to  a  resonant  acoustic  sound  field.  It  was  believed 
that  the  data  obtained  from  the  flat  plate  system  would  make  more  and  simpler 
experimental  facts  available  than  had  previously  been  obtained  with  the  hori¬ 
zontal  tube  heat  transfer  apparatus  at  Georgia  Tech.  This  reasoning  was  par¬ 
ticularly  valid  if  natural  convection  was  complicating  the  experimental  measure¬ 
ments.  Unfortunately,  experimental  difficulties  were  encountered  with  the  flat 
plate  system  and  it  was  impossible  to  obtain  useful  information  from  it.  Some 
of  the  reasons  for  the  experimental  problems  have  been  given  in  Chapter  III. 
Preliminary  visual  exploratory  studies  with  a  tube  geometry  indicated  that  free 
convection  was  not  complicating  the  secondary  flow  and,  therefore,  studies  were 
carried  out  with  this  system. 

It  is  interesting  to  note  that  the  theoretical  flat  plate  analysis,  limited 
for  the  present  time  to  flows  without  heat  transfer,  gives  some  very  interesting 
and  useful  results  and  even  indicates  some  plausible  explanations  of  some  of  the 
previous  experimental  data  on  tube  flow.  The  analytical  study  also  indicates 
how  a  similar  treatment  of  the  circular  geometry  can  be  attempted. 

The  flow  field  shown  in  Figure  7  presents  a  plausible  explanation  of  the 
data  shown  in  Figures  1  and  2  and  is  currently  being  studied.  The  maximum  local 
heat  transfer  rates  occur  where  the  main  vortex  and  the  main  stream  impinge  upon 
the  wall,  i.e.  at  the  velocity  loops.  The  minimum  rates  occur  where  both  the 
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vortex  and  the  main  stream  separate  from  the  wall,  i.e.  at  the  velocity  nodes. 

Another  interesting  outgrowth  of  this  study  was  the  correlation  of  the 
threshold  data  of  Eastwood  [2].  Eastwood  measured  the  local  heat  transfer 
rate  for  forced  convection  in  a  horizontal  isothermal  tube  with  an  impressed 
resonant  acoustic  field.  He  defined  the  sound  threshold  level  to  be  the  sound 
level  below  which  no  noticeable  change  in  the  local  heat  transfer  rate  occurred. 
If  this  is  interpreted  in  terms  of  the  flow  field  which  exists,  then  it  should 
follow  that  a  certain  magnitude  of  disturbance  in  the  flow  would  be  necessary 
to  produce  an  observable  change  in  the  convection  heat  transfer  characteristics 
of  the  tube.  From  the  present  theoretical  analysis  the  magnitude  of  this  dis¬ 
turbance  should  be  characterized  by  the  parameter  m/m^.  Eastwood's  threshold 
data  were  correlated  at  a  value  of  equal  to  63  for  the  range  of  average 

i4ach  number  between  0.005  and  0.017.  For  Mach  numbers  less  than  0.005  the 
effect  of  free  convection  became  appreciable  and,  therefore,  below  this  value 

was  not  a  sufficient  description  of  the  f)ow. 

A  point  which  should  be  emphasized  again  is  that  the  analytical  solution 
presented  herein  for  channel  flow  is  only  valid  for  values  of  less  than 
approximately  10.  However,  the  fact  that  Eastwood's  threshold  data  correlated 

for  a  value  of  of  63  indicates  that  the  restriction  may  be  too  severe. 

0 

It  is  interesting  to  note  that  if  the  through  flow  is  reduced  to  zero  then 
the  time-mean  solution  is  the  same  as  that  predicted  by  Westervelt  [6].  The 
other  extreme  case,  that  of  no  sound,  reduces  the  solution  to  one  for  fully 
developed  laminar  channel  flow.  This  condition,  however,  is  not  required  by  the 
analysis  presented  here  and  should  be  considered  as  a  curin<«lty  rather  than  a 
necessity. 
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Preliminary  results  of  the  quantitative  visualization  experiments  with  tube 
flow  are  shown  in  Figure  15.  The  comparison  with  the  theoretical  solution  for 
channel  flow  is  interesting.  Also  the  lack  of  scatter  in  the  experimental  data 
indicates  that  the  parameter,  WUT,  is  apparently  a  correlating  parameter  for 
tube  flow  as  well  as  for  channel  flow.  As  mentioned  previously  this  parameter 
correlated  Eastwood's  heat  transfer  data  for  tube  flow. 

Although  differences  between  the  flow  patterns  of  air-oil  vapor  and  air- 
tobacco  smoke  were  not  apparent,  qualitative  differences  did  exist.  Since  the 
particle  size  of  tobacco  smoke  is  much  smaller  than  that  cf  vaporized  oil,  only 
tobacco  smoke  was  used  to  obtain  the  quantitative  data.  Typical  flow  patterns 
are  shown  photographically  in  Figure  10. 

It  is  anticipated  that  the  solution  for  the  circular  duct  will  be  completed 
by  the  first  of  the  year  and  that  further  improvements  in  the  experimental 
techniques  will  permit  an  accurate  verification  of  this  solution. 
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Figure  1.  Local  Nusselt  Numbers  Versus  x/L  for  Various  Sound 
Pressure  Levels  and  221  cps. 
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Figure  3*  Local  Nusselt  Numbers  Versus  x/L  for  Various  Sound 
Pressure  Levels  and  222  cps. 
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PUgure  5*  Time  Dependent  Velocity  Profiles. 
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Figure  6.  Time  Mean  Velocity  Profiles. 
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Ficure  9.  Sound  Pressure  Level  Versus  Mean  Through  Flow  Mach  Number 
for  Various  Vortex  Sizes. 


FLOW 

(a)  1440  cp* 

M  -  1  X  10"* 


(b)  1600  cpt 
M  -  1  X  10"* 


(c)  1120  cpt 

M-O.Sx  10”* 


Figure  10.  Qualitative  Smoke  Patterns  from  Preliminary  Investigations 
for  Tube  Flow. 
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FLOW 

FREQUENCY  .  1000  ep* 
MACH  NUMBER  -  0.0001 


Figure  11. 


Acoustic  Vortices  for  a  Rssor.ar.o  Recoar-rular  D-uco. 
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I'ij'utv  It?.  Ucticinutlc  Dlu('mm  oI‘  Flow  VlsuuJ  Izutlon  Equipmcn't. 


Photcgraph  of  KLow  ^sualizaticn 
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Figure  15-  Canparison  of  Experimental  7S  for  Tube  with  Analytical  7S  for 
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